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Abstract 

We study the time evolution of wave packets of noninteracting electrons in a two- 
dimensional periodic system in the presence of magnetic and electric fields. The 
model includes consistently the coupling between Landau levels as well as the peri- 
odic and external field contributions. It is shown that the electric field produces lo- 
calization along its longitudinal direction, this effect is a physical consequence of the 
quantization of the energy levels with respect to the longitudinal quasi-momentum. 
For incommensurate orientations of E relative to the lattice the wave packet be- 
comes localized in both direction. 
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Localization of waves in quasiperiodic system has been the subject of great in- 
terest in condensed matter and related areas in physics. An electron subject to 
a sum of potentials which are periodic but incommensurate with one another 
represents that is intermediate between a random and a perfectly or- 

dered system. In an incommensurate situation there can be both localized and 
extended states, with the intermediate possibility of a metal-insulator transi- 
tion as the energy and-or the strength of the potential are varied. Aubry and 
Andre [1] presented a one-dimensional model that exhibits this behavior; the 
dynamics of their model is governed by the Harper [2]. The same model also 
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arises in the the study of an electron moving in a two-dimensional lattice sub- 
ject to a magnetic field, both in the limiting cases of either a very large [2] or a 
very small [3] potential strength compared to the magnetic field intensity. The 
problem of electrons moving under the simultaneous influence of a periodic 
potential and a magnetic field has been discussed by many authors [2,3,4,5]; 
the spectrum displays an amazing complexity including various kinds of scal- 
ing and a Cantor set structure [6]. The addition of an electric field leads to 
interesting new phenomena that makes its analysis worthwhile. For example, 
as the strength of the electric field increases the longitudinal quasi-momentum 
is quantized leading to the appearance of a "magnetic Stark ladder" in which 
the bands are replaced by a series of quasi discreet levels [7,8]. The object of 
the present work is to study the time evolution of electron wave packets in 
a periodic two dimensional system subject to external electric and magnetic 
fields. Previous studies of this problem were usually restricted to a single band 
limit derived from a tight binding and/or the Peirles approximations [9,10]. 
Here we presents a model which consistently includes the coupling between 
Landau levels as well as the periodic and external field contributions. This 
treatment becomes essential if one would like to analyze physical realistic sit- 
uations. The study of these phenomena has become experimentally accessible 
with the recent developments in the fabrication of antidot arrays in lateral 
superlattices by ion beam and atomic force lithography [11,12]. For example, 
the observation of the Hofstadter Butterfly spectrum through the measure- 
ment of the magnetoresistance and Hall conductance in artificial arrays of 
anti quantum dots has just been achieved [13]. 

We consider the motion of an electron in a two-dimensional periodic po- 
tential V (x, y) subject to a uniform magnetic field B perpendicular to the 
plane and to a constant electric field E, lying on the plane according to 
E = E {cosQ, sin9), with 9 the angle between E and the lattice x— axis. The 
dynamic is governed by the Schrodinger equation 



the vector and scalar potentials are selected to yield: V x A = Bk, and 
E = — V0 — dA/dt. For simplicity we shall consider a periodic square lattice: 
V (x, y) = Uq [cos {2irx/a) + A cos (2iry/a)} . Two characteristic parameters of 
the system, the cyclotron frequency uo c = eB/m and the magnetic length 



Ib = \ffr/( e B) set the reference scales for energy and lengths respectively. 
Following the formalism developed in references [7,8], we introduce a canonical 
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transformation to new variables Q, P, acoording to 
Qo — t , 

eQi = P y + eA y + mE x /B , 
eQ 2 =Px + eA x - E x t , 

where = + r • E, A x = A x + By + E x t and A y = A y + Bx + E y t. It 
is easily verified that the transformation is indeed canonical, the new vari- 
ables obey the commutation rules: [Q 0) -Po] — —ifi, and [Qi,-Pi] = [C&j-f^] = 
ihB; all other commutators being zero. The inverse transformation gives x = 
1% (eQi — P2) /h and y = l\ (eQ2 — Pi) /H. The operators (Pq, Q2, P2) can be 
identified with the generators of the electric-magnetric translation simetries 
[14,7]. Final results are independent of the selected gauge. From the operators 
in Eq. (2) we construct two pairs of harmonic oscillator-like ladder operators: 
(ai,a\), and (02,02) with : 

obeying: [ai,a|] = [02,02] = 1, and [01,02] = [01,02] = 0. Using these vari- 
ables, the Schrodinger equation can be recast as P IV') = H \ip) with the 
Hamiltonian given by 

H = hu c (^ + a i a i) +y~el B (Ea 2 + E*a 2 ) , (4) 

here E = (E x + iE y )/y/2, and V(x,y) = V(Q X - P 2 ,Q 2 - Pi) has to be 
expressed in terms of the ladder operators using Eq. (3). The properties of the 
system depends on five dimensionless parameters: (i) A modulates the relative 
potential amplitudes in x and y directions, A > 1 (A < 1) corresponds to a 
weaker modulation along the x (y) direction, (ii) a^ 1 = <p/(po = a 2 B/(h/e) 
is the number of flux quanta in a unit cell of area a 2 , (iii) 9 is the relative 
angle between the electric field direction and the lattice The two other 

parameters K and p are the potential strength in units of Hou c , and the electric 
field in units of Uo/(ea), hence: (iv) K = Uo/hu c , and (v) p = ea\E\/Uo, 
respectively. As an appropriate set of basis functions we choose the eigenstates 
of the number operators ajai and a\a2 

a\ai \p,, v) = p \p, v) , a\a 2 \p,v) = v\p,v) , (5) 

where p labels the Landau levels and \p, v) — \p) <g> \ v). These states are also 
eigenvectors of the angular momentum operator J = (P 2 + Q\ — P% — Ql)/2: 
J \p, v) = h(p — u) \p, u) . In the coordinate representation the wave function 



P = t hd t + e4>+-(\E\/B) 2 , 
Pi =p x + A x -mE y /B, 

P 2 =Py + eAy - Eyt, (2) 
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for this state takes the form 
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(6) 



where L v ~^ are the generalized Laguerre polynomials, the previous formula 
is valid for /i < v, if /i > v the indices fi and v must be exchanged on 
the right hand side. These wave functions are centered in the origin and its 
mean square radius is given by (fi, u\r 2 \ fi,u) =81% (fi + v + 1) . The initial 
form of the wave packet has been selected as the lowest Landau level with 
zero angular momentum: ^0,0 ( x ,y)- The evolution of the system is obtained 
applying the operator O (t) = exp (—itH) to the initial state i.e. \1/ (t,x,y) = 
{x,y \0 (t) \ 0,0) . The Hamiltonian in (4) is represented by a matrix EI with 
elements evaluated in the base given in Eq. (5), the matrix elements can be 
worked out as 



(n',i>'\H\n, v) 



+ 



Un 



(7) 



where D^^ (z) = (ji' e zai z * a i fi} = Vy,^ (\^2Rez, y/^lmz^j are the matrix 

elements of the coherent state operator e ( za i~ z * a i) ; they can be expressed using 
the same functions that determine the wave function profiles (Eq. 6). The 
Hamiltonian in (7) is diagonalized numerically by the unitary matrices U that 
are constructed from the column eigenvectors of H. In this base the evolution 
operator takes the form O (£) = U 1 " exp (— itE)V, where E is the diagonal 
matrix containing the energy eigenvalues. The time evolution of the system is 
given by 

00 

*(t,x,y)= £ C^(t)^ u (x,y), (8) 

H,u=0 

where C (t) = O(t)C(0), with the initial state C£ (0) = 5 M ,o<W A global 
characterization of the dynamic evolution of a wave packet is provided by 
displaying 3 D graphics of the probability density p (t, x,y) = \ty (t, x, y)\ 2 , as 
a function of the space-time variables. The localization properties of the system 
are characterized by the evolution of the mean square displacement u 2 (t). 
Both the electric field and the A parameter can induce anisotropy transport 
properties, so its is sometimes convenient to separately refer to the x and y 
variance: a\ (t) = (x 2 } t - (x) 2 , and a 2 (t) = (y 2 } t - (y) 2 , with a 2 , = a 2 x + a 2 . 
As an additional criteria we consider the Shannon information entropy [15], 
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t = Or 



(a) E = 



t = 2t 



t = Or 



(b) E = E i 

t = 4r 




Fig. 1. Time evolution of an initially localized wave packet in a 2D periodic asym- 
metric potential (A = 2; lower modulation amplitude in the x direction) for: (a) 
E = 0, and (b) an electric field pointing in the x direction, with an intensity 
p = ea\E\/Uo = 0.1. Time is measured in units of r = 2it/lu c . Here and throughout 
the paper the value of the other parameters are a = 1/2, K = 5.5. 



defined by 



S(t) 



p(t, r) In p(t, r) d 2 r. 



(9) 



This quantity characterize the spreading of the probability density distribu- 
tion, it provides a measure of the spatial derealization of the wave packet. 
The lower this quantity is, the more localized the electron is. 

In what follows we analyze the propagation of wave packets and its properties 
for some specific examples. Using the Peirls tight-binding model Nazareno 
and Brito [9] studied the electron wave propagation, they found that for a 
rational the electron propagates ballistically, whereas for irrational a the wave 
remains localized. Using our model that incorporates inter-Landau couplings, 
we have verified that these results remain valid; details are not presented due 
to space limitations. In what follows we shall restrict to the case of a — 1/2 
and fix the parameter K to K = 5.5. The Hamiltonian in (7) is diagonalized 
numerically considering a base of p = 4 landau levels and v = 850 excitations, 
this selection proves to give good convergence. This was verified by checking 
that the participation of states with higher values of p and v remains smaller 
than one part in 10 7 . 

Directional localization induced by the lattice anisotropy and the 
electric field. If the lattice potential is anisotropic (A 7^ 1) electron trans- 
port is expected along the direction of weaker modulation. Fig. la presents 
the time evolution of a wave packet initially localized, in the absence of elec- 
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Fig. 2. (a) Time evolution of the a 2 and a 2 variances for the wave packet rep- 
resented in Fig. lb, the wave packet is dynamically localized in the longitudinal 
direction and spreads ballistically in the direction of stronger modulation a 2 ~ t 2 . 

(b) Energy subspectra as a function of the parameter p (p = ea\E\/Uo). The trans- 
verse quasimomentum is fixed to k y = and the longitudinal quasimomentum takes 
all possible values in the first Brillouin zone, k x £ (—ir/2a,+ir/2a). The first two 
Landau levels are shown, a = 1/2, so every Landau level splits in two bands. A 
transition from extended to localized states is observed as the electric field intensity 
is increased. 

trie fields, and for a lower modulation along the x-direction (A = 2). Clearly 
as expected, ballistic transport takes place along the direction of lower mod- 
ulation. If an electric field in the y direction is applied, the characteristic 
properties of ballistic transport along the x-axis would remain valid. Instead, 
let us consider that an electric field pointing in the direction of lower mod- 
ulation (rr-axis) is switched-on, for sufficiently strong field a transition takes 
place and the electrons are re-directed into the direction of stronger modula- 
tion and simultaneously localized in the direction of lower modulation (x-axis). 
This behavior is displayed in Fig. lb, we observe that initially the effect of the 
lower potential modulation dominates and the wave package spreads in this 
direction until it reaches a maximum size (at a time t ~ 10 r, r = 2ty/uj c ), 
after this time the diffusions along the longitudinal directions stops, with the 
mean square displacement a 2 showing strong oscillations; simultaneously the 
package spreads along the transverse direction (y axis). This behavior is cor- 
roborated by the behavior of the x and y variance. In Fig. 2a. a 2 and a 2 are 
plotted as function of time, the evolution in the transverse higher modulation 
direction is ballistic, whereas localization in the longitudinal direction with 
a strong oscillatory behavior for a 2 is observed. Hence the system has un- 
dergone a metal-insulator transition in the x (longitudinal) direction and a 
insulator-metal transition in the y (transverse) direction. The high amplitude 
of the a 2 oscillations have its origin in the competing contributions between 
the modulation of the periodic potential and the electric field. The change 



6 



in the direction of transport, forcing the electron to overcome the stronger 
potential modulation, is related to the behavior of the spectrum and eigen- 
functions of the Hamiltonian (7). Fig. 2b presents the energy subspectra as 
a function of the electric field strength for the first two Landau levels, no- 
tice that for a = 1/2 every Landau level splits in two bands. To generate 
the subespectra, the transverse quasimomentum is fixed to a value k y = 0, 
while the longitudinal quasimomentum takes all possible values in the first 
Brillouin zone, k x G (— 7r/2a, +7r/2a). It is observed a transition from wide to 
extremely narrow minibands, simultaneously the eigenfunctions change from 
extended to localized with respect to the longitudinal quasi-momentum. This 
extremely narrow band corresponds to the generation of quasi discreet levels 
or "magnetic Stark ladder" [7,8], and leads to the longitudinal localization. 

A phenomena similar to the one presented here, has been reported by Ket- 
zmerick, Kruse and Springsguth [16] . In that case, the transition that induces 
transport along the direction of stronger modulation and localization in the 
direction of weaker modulation is related to avoided band crossing that pro- 
duces transitions from band to discreet levels or vice versa as the strength of 
the periodic potential is increased. In our present example the transition is 
related to the change from wide to extremely narrow minibands induced by 
the electric field. 

Commensurate vs. incommensurate directions of the electric field Fi- 
nally, we consider the effects produced by the commensurability of the electric 
field direction. We refer to a commensurate direction if the electric field is ori- 
ented in such that a way that the angle fulfills the condition tan 9 = 1712/1711, 
with uii and m<i relatively prime integers, this condition ensures that a spatial 
periodicity is preserved along the longitudinal and transverse direction of E. 
In the rest of the paper we shall restrict to a symmetric lattice (A = 1). (i) 
Let us first consider the simplest commensurate case 9 = 0. The wave package 
spreads along the transverse direction, Fig. 3a. The mean square a x presents 
small amplitude oscillations, while a 2 , shows a ballistic growing, Fig. 4a. (ii) 
Incommensurate orientation: 9 = tan -1 aa- Fig. 3b clearly shows that the 
wave packet is dynamically localized in both directions, within a finite re- 
gion of the lattice. In Fig. 4b it is shown that both a 2 and a 2 oscillate. It 
is interestig to compare the periods of these oscillations: the period of a 2 is 
T x m 6.9r and a 2 has a period T y 11. It. The ratio of this two periods is 
similar to the golden ratio T x /T y ~ E y /E x = (\/5 — l)/2, which seems to be a 
consequence of the relation between the two components of the drift velocity 
v D — E (sin 9, — cos 9) /B. Along the x direction the electron will take a time 
T x = qa/vc sin 9 to reach the extreme of the unit cell, while on the y axis 
we have T y = a/vp sin 9, giving the correct relation between the two periods. 
Each time the electron reaches the border of a unit cell it reflects showing 
an oscillation in the mean square dispersions. The oscillations are damped 
because of tunneling. 
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(a) = 
t = Or t = 3r 



(b) = tan" 1 a G 
t = Or t 'dr 




Fig. 3. (a) Wave packet evolution when the electric field is oriented in an com- 
mensurate direction 9 = 0, the electron is localized in the longitudinal direction, 
(b) Wave packet evolution when the electric field is oriented in an incommensurate 
direction 9 = tan -1 olq, = [(\/5 — l)/2]. The wave packet is localized in both 
directions. In both examples A = 1, and p = 0.1 

The comparison of the commensurate (Fig. 3a) with incommensurate (Fig. 3b) 
cases clearly shows the very different transport properties of the electron de- 
pending on the relative orientation of the electric field with respect to the 
lattice. The behavior of the Shanon entropy is consistent with these results 
(Fig. 4c), we observe that for small times the evolution of S is common for both 
cases, however for t > 3r the growth of S(t) corresponding to the irrational 
orientation is considerable reduced as compared to the rational case. 

The localization for incommensurate directions of the electric field can be un- 
derstood with the following argument. For the electric field orientation tan 9 = 
m 2 /mi, the periodicity is preserved along the longitudinal and transverse di- 
rection but with an extended lattice of side dimensions b = a \jm\ + m\. If 
additionally, the number of flux quanta in each of the new unit cells is given 
by a rational number, i.e. o~ x = b 2 B/(h/e) = p/q, then the electric-magnetic 
translation symmetries are preserved in the extended lattice of dimensions 
qb x b = q (m 2 +m|) a 2 [7,8]. With the new conditions, the on-site energies of 
the original lattice are no longer the same, in order to coherently propagate 
the electron has to tunnel a distance proportional to the new lattice dimen- 
sions q (m\ + m|) a 2 . Clearly as m\ or m 2 increases the wave packet diffusion 
will be inhibited. For irrational orientations of the electric field, we can con- 
sider the rational approximant tan 6 = /w!^\ then the system becomes 
quasi-periodic and the electrons are not allowed to propagate coherently thus 



8 



8 = 

= tan -1 a 



5 10 15 
t 

Fig. 4. (a) Time evolution of the and variances when the electric field is 
aligned in a commensurate direction 6 = 0. (b) Time evolution of cr^ and o~y for an 
incommensurate direction 6 = tan -1 ag of E. In both examples A = 1, and p = 0.1. 
For the rational orientation of E the electron is only localized in the longitudinal 
direction, whereas for 6 = tan -1 ac it is localized in both directions, (c) Time 
evolution for the Shannon entropy for the two orientations of E. 

leading to localization. 

We have analyzed the time evolution of electron wave packets in a two- 
dimensional periodic system in the presence of magnetic and electric fields. 
The dynamics is governed by the effective Hamiltonian in (7), which consis- 
tently includes the coupling between Landau levels as well as the periodic and 
external field contributions. The inclusion of an electric field (for commensu- 
rate orientations) yields localization in the longitudinal direction; propagation 
takes place along the transverse direction. This effect is a physical consequence 
of the quantization of the energy levels with respect to the longitudinal quasi- 
momentum. Transition from extended to localized states are observed as tan 9 
changes from rational to irrational values. It is expected that these phenom- 
ena should be experimentally accessible using lateral superlattices on semi- 
conductor heterostructures, for example by measuring assymetry effect on the 
magnetoresistance and Hall conductance as the electric field orientation with 
respect to the lattice axis is varied. 
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